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Abstract 

The large n behaviour of the hypergeometric polynomial 

3^2 I i i , " 

\i - n, i - n 

is considered by using integral representations of this polynomial. This 
3F2 polynomial is associated with the Catalan-Larcombe-French sequence. 
Several other representations are mentioned, with references to the liter- 
ature, and another asymptotic method is described by using a generating 
function of the sequence. The results are similar to those obtained by 
Clark (2004) who used a binomial sum for obtaining an asymptotic ex- 
pansion. 
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1 The problem 

Find the large n asymptotics of 

/<"> -^(jIm .„;-») <"> 

Peter Larcombe conjectured that linin^oo f(n) — 2 and Tom Koornwinder gave 
a proof, based on dominated convergence. See for details of the proof |S] , where 
also a different representation of f(n) is considered in the form 



7/i 



/(n)=2» 3 F 2 ( n *' ; I )• (1.2. 



5 ~ 1 
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The equivalence of these two forms follows from a quadratic transformation of 
the 3-F2— functions as given in [5J Ex. 4(iv), p. 97], that is, 

a, b, c 

■■J--> I ; z 

1 + a-b, 1 + c-c 



(1 - z)- a 3 F 2 



(1.3) 

\a, \ + ^a, 1 + a — b — c — 4z 



1 + a — b, 1 + a — c 



with a = — n, 6 = c = i, and z = — 1. Another form is given by (see Eq. 
(A2)]) 

-n, — n, 
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In 0] an asymptotic expansion of \f{n) has been derived. The asymptotic 
analysis is based on the representation 

P ,4 E (1 . 5) 

p-\-q—n y ' 

By using the relation 



(2n)! = 2 2 "n!(±)„, ? i = 0,l,2,..., (1.6) 

it is straightforward to verify that (|1.5|) can be written as 

p _ 2 4 " (2 )p(2")p(2")"~p( 2^ n ~P (17) 
™ n! -^-^ p! (n — p)l 



By using 



it follows that 



(a) w _ fe = (-l) fc n (a) " , (1.8) 
(1 - a - n)k 



that is, 



p _ 2 4 "( 2 )„( 2 ) ra ( "-)p( 2 )p(2)p Q\ 

"~ n\n\ ^ !) p| (J -„) p (i -n) p ' (L9) 



n) n\ \ -k — n. -h — n / v ' 



nln\ \ i - n, ^ — n 

which gives the relation with /(n) by using l|l.l|l 



Pn = ^%^L/(n) = ( 2n Xf(n). (1.11) 

~t.\ n.\ 
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The numbers P n are for n — 0, 1, 2, , . . . known as the elements of the se- 
quence (A053175) {1, 8, 80, 896, 10816, . . .}, called the Catalan- Larcombe- French 
sequence, which is originally discussed by Catalan See the On-Line Ency- 
clopedia of Integer Sequences http:/ /www. research. att.com~njas/sequences/ . 

In this paper we derive a complete asymptotic expansion of the numbers 
P n by using integral representations of the corresponding 3 F2— functions. Our 
results are the same as those obtained by Clark 0] , who used the binomial sum 
in (|1.5|) without reference to the 3F2 — functions. 

2 Transformations 

We derive an integral representation of the 3F2 — function of (|l.l|l by using several 
transformations for special functions. We start with the beta integral 




(2.12) 



and use it in the form 




0,1,..., n. (2.13) 



We substitute this in the representation of the 3F2— function in i|l.l|l 




(2.14) 



This gives after performing the k— summation 




(2.15) 



We substitute t = sin 2 (0/2) and obtain 




(2.16) 



We apply a quadratic transformation (see ^ Eq. 15.3.26]) to obtain 




(2.17) 



and use the representation of the Legendre polynomial 




(2.18) 
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Figure 1: The domain of integration of the integral in 12,211 and subdomains A, B, 
C and D. 



This follows from [3 Eq. (22.3.8)] and gives 



f W = m m / sin " 9 p « 77 ) dd - ( 2 - 19 ) 



Next, consider (see jS] p. 204]) 

P„(z) = i^ (z + ^ 2 -lcosV>)" #, n = 0,1,2,..., (2.20) 



which gives the double integral 
n\ n! 

71-2 (i)" (s)n 7o 7o 



„. . n!n! r / 1 + cos 6 cos ji \ , , 
/ B ) = -J7TwTr / S ( 2 - 21 ) 



3 Asymptotic analysis 

The landscape of the integrand in l|2.21|l shows peaks at the boundary points 
(0,0) and (71", 71"), where it assumes the value 1. Along the interior lines 9 — hir 
and tp = |tt the integrand has the value 2~™. Inside the squares A and C, see 
Figure^ the value of the integrand is between 2~™ and 1, in the squares B and 
D it is between and 2~ n . In addition, the contributions from A and C are the 
same, and also those from B and D are the same. 

From an asymptotic point of view it follows that the integral over the full 
square equals twice the integral over A, with an error that is of order C(2 _n ), 
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while the total integral is of order 0(1), as n is large. Hence, we concentrate on 
the integral over A, and write for large values of n 



„, . „ 2- n n\n\ 
fin) = 2- 



71-2 (i)« (3)7 



7T /■ tt7T 



( 1 + cos 9 cos ip) n dO dijj + E n 



(3.22) 



where E n = 0(2 "). Next, we neglect E n and put u = sin(9/2),v = sm(ijj/2), 
and obtain 

^ (i)n(i)™Jo Jo Vl-u 2 Vl-w 2 

(3.23) 

For the integrals in i|3.22|l and l|3.23|l asymptotic expansions can be obtained 
by using Laplace's method for double integrals; see ^3 § VIII. 10]). In our case 
a simpler approach is based on neglecting a part of square A by introducing 
polar coordinates 

u = rcost, v = rsint, < r < ±V2, < t < ±tt. (3.24) 
This gives (again we make an error in the integral that is of order 0(2 - ™)) 
8n!n! f n/2 f ^ (1 - r 2 + 2r 4 cos 2 t sin 2 t) n rdrdt 

A n ) 2n\ m / / — / ( 3 ' 25 ) 

n(2)n( 2 )nJo Jo ^(1 - r 2 cos 2 t) (1 - r 2 sin 2 i) 

We change r 2 into r, and obtain 

An\n\ r /2 fi (1 -r + 2r 2 cos 2 < sin 2 £)" dr (ft 

/W 2m m / / — / ( 3 - 26 ) 

n y^^hJo Jo y'(l-rcos 2 i)(l-rsin 2 i) 

First the standard method for obtaining asymptotic expansions of a Laplace- 
type integral can be used (for the r— integral). The second step is done by 
integrating the coefficients of this expansion with respect to t. 

For the r— integral we transform the variable of integration by putting 

w = - In (l - r + 2r 2 cos 2 t sin 2 t) . (3.27) 

This mapping is one-to-one for re [0, 5], uniformly with respect to t S [0, ^tt], 
with corresponding w— interval [0,U»o], where wq = w(h). 
We obtain 

/(")- 2fl ' 7T-, / / e~ nw F(w,t)dwdt, (3.28) 

n \2)n {2~)n Jq Jo 

where 

F(w, t) = 1 (3.29) 

(l-rcos 2 i)(l-rsin 2 t) dw 
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4 Asymptotic expansion 

We obtain the asymptotic expansion of w— integral in (|3.28|) by using Watson's 
lemma (see § 1.5]). 

The function F(w,t) is analytic in a neighborhood of the origin of the 
w— plane. We expand 

oo 

F(w,t) =^ Cfe (i)u, fe (4.30) 

fc=0 



where 



and substitute this expansion in (|3.28l) . Interchanging the order of summa- 
tion and integration, and replacing the interval of the w— integrals by [0, oo) (a 
standard procedure in asymptotics) we obtain 

i> i v An\n\ ^ fc! . 

C fe = / Cfc(t)dt, fc = 0,1,2,... . (4.32) 
Jo 

The coefficients Ck(t) can be obtained by the following method. First we need 
the inverse of the transformation defined in 13.27fl . That is, we need coefficients 
bk in the expansion 

oo 

r(w)=J2 bk ^ wk - ( 433 ) 

We can find r(w) from (|3.27() as a solution of a quadratic equation, with the 
condition r(w) ~ w as u> — > 0, that is, &o(i) = !■ However, we can also 
differentiate (|3.27|l with respect to r and substitute the expansion l|4.33|l . and 
solve for the coefficients bk(t). When we have these coefficients we can expand 
F(w,t) of and find c k {t). 

The first few coefficients Ck(t) are 

co(t) = 1, 

ci(i) = ±(-l + 8s 2 -8s 4 ), 

c 2 (i) = |(1 - 28s 2 + 220s 4 - 384s 6 + 192s 8 ), 

c 3 (i) = + 92s 2 - 1628s 4 + 10752s 6 - 24576s 8 + 23040s 10 - 7680s 12 ), 

c 4 (i) = 3§j (1 - 280s 2 + 10024s 4 - 130848s 6 + 773904s 8 - 2054400s 10 + 

2691840s 12 - 1720320s 14 + 430080s 16 ), 
c 5 (t) = 3^q(-1 + 848s 2 - 55328s 4 + 1259040s 6 - 13396560s 8 + 

73983360s 10 - 215329920s 12 + 349224960s 14 - 

319549440s 16 + 154828800s 18 - 30965760s 20 ) 

(4.34) 
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where s = sin 2 1. For the corresponding Ck we have 

Co = |tT, C\ = 0, C 2 = g7T, C 3 = |7T, C4 = ^7T, C 5 = gj7T. (4.35) 

As a next step we can replace in (|4.31|l the ratios n!/(^)„ by the asymptotic 
expansion 

nl T(n+1) 7fc 

— 7r — - 7- r-j V7rn x 



-r— - \ r- - x -'T-T. ( 4 -36) 

where 



^° 1' 71 g j 72 128 , 73 1024 , 74 32768' ^ 5 262144' 

(4.37) 

This finally gives 

„. , / 1 17 207 14875 352375 \ on , 



5 An alternative method 



The numbers P„ were proposed as "Catalan" numbers by an associate of Cata- 
lan. They appear as coefficients in the series expansion of an elliptic integral of 
the first kind 

K(k)= 1 dt, (5.39) 

Jo VI - k 2 sin 2 t 

which is transformed and written as a power series in k (through an intermediate 
variable); this gives a generating function for the sequence {P n }- For details we 
refer to [S]. 

In [5] a generating function for the numbers P n is given in terms of the square 
of a modified Bessel function, and we use this approach to obtain an asymptotic 
expansion of fin). See also |7| for details on this generating function. 

We consider numbers F n defined as coefficients in the generating function 



e w / 2 I (w/2)] =J2 F " 1 



(5.40) 



n=0 



By considering the relation of the Bessel function with the confluent hypergeo- 
metric functions (see P Eq. 13.6.3]), 



e z !Az) 



far 



2 ■ 2z 



T{u + 1) 

we can write (|5.40|l in the form (see also Eq. 13.1.27]), 





to 


iFi (j ; w^j 


= e 2w 



1*1 



J2F n w n . 



n=0 



(5.41) 



(5.42) 
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This gives the representation for F n : 

{o)k (o)n.— 



fc=0 



k\k\ (n-fc)!(n- fe)! 



By using (|1.8|l it follows that 



F„ 



).! 77.1 



EC- 1 )- 



; (-rc)fc(-re)fc(^)fc 

{\-n) k k\k\ 



or 



(2)" p 



It follows from (|T^|) that 



From H5.40JI we obtain 



fin) 



-71, — n, 5 
l,|-n 



-1 



? ! n! «.! 



2" (2)" 



1 

2ni 



„,n+l 



cz/u; 



1 



2ni Jr w 



—r h(w) dw, 
1+1 v ' 1 



where 



h(w) 



- w/2 I (w/2) 



iF 



-w 



(5.43) 

(5.44) 
(5.45) 

(5.46) 

(5.47) 
(5.48) 



and the contour C is a circle around the origin, or any contour that can be 
obtained from this circle by using Cauchy's theorem. The main contribution 
comes from the saddle point of , that is from w — Wq = n/2. 

In the standard saddle point method (see ^3 § II. 4]) a quadratic transfor- 
mation is used to bring the main part of the integrand in the form of a Gaussian. 
We can obtain the same expansion by just expanding the function h(w) (which 
is slowly varying for w > 0) at the saddle point. 

First we expand (see P Eq. 13.4.9]) 



/.'.i. 1 j ; -w 



and next 



= y^kjw - wa) k 

fc=0 



(-in*)* 



k\k\ 



2/ ~iF 



\ + k 
1 + fc ' 



-w 



h(w) =J2 A k( w ~ w o) k - 



(5.49) 
(5.50) 



We substitute this expansion in the second integral in (|5.47|) and obtain the 
convergent expansion 



~ l f e 2w 



(5.51) 
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The functions <E>/c can be evaluated by using the recursion formula (which easily 
follows from integrating by parts) 

2 n 

* fc = -i(fc-l)($ fc _ 1 +i« $ fc _ 2 ), $ = ^, $i=0. (5.52) 

An asymptotic expansion can be obtained by using a well-known expansion 
for ak defined in H5.49fl . We have (as follows from ^ Eq. 13.5.1]) 

lFl ( a ; -x] ~ x-^L f {a)m(1 + a ~ c ^ , x ^ +00 , (5.53) 

v 7 K ' m=0 

from which we can obtain expansions for <Xk and Afc for large values of wq = n/2. 
By using these expansions in (|5.51|) we obtain an expansion for F n , and finally 
for f(n by using (|5.46l) . This expansion is the same as the one in l|4.38() . 
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